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The Discrete Analog of the Malgrange-Ehrenpreis Theorem 

Doron ZEILBERGER^ 

In fond meomory of Leon Ehrenpreis 

One of the landmarks of the modern theory of partial differential equations is the Malgrange- 
Ehrenpreis [E][M] theorem (see [Wi]) that states that every non-zero linear partial differential 
operator with constant coefficients has a Green's function (alias fundamental solution). Recently 

Wagner [W] gave an elegant constructive proof. 

In this short note I will state the discrete analog, and give two proofs. The first one is Ehrenpreis- 
style, using duality, and the second one is constructive, using formal Laurent series. 



^^ ■ Let Z be the set of integers, and n a positive integer. Consider functions f{mi, . . . ,mn) from Z"" 

to the complex numbers (or any field). A linear partial difference operator w^ith constant 
^^ . coefficients V is anything of the form 



Vf{mi,... ,m„) := ^ Caf{mi + ai,..., rUn + a„) 



where A is a finite subset of Z" and a = (oi, . . . , a„), and the c^ are constants. 
^ For example, the discrete Laplace operator in two dimensions: 

o. ^ 

°^ ; /(mi, ma) -^ /(wii, mg) - -(/(m-i + 1, "^2) + /("^i - 1, "12) + f{mi,m2 + 1) + f{mi,m2 - 1)) • 

^/ 

(^ ' The symbol of the operator V is the Laurent polynomial 



-r[Zi, . . . , Zn) — / ^ CaZ^ • • • , Z„" 

The discrete delta function is defined in the obvious way 

5(mi,...,m„) = 1^' ^Hmi,...,mn) = (0,0, ...,0); 
[0, otherwise. 

Note that the beauty of the discrete world is that the delta function is a genuine function, not a 
"generalized" one, and one does not need the intimidating edifice of Schwartzian distributions 
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We are now ready to state the 

Discrete Malgrange-Ehrenpreis Theorem: Let V be any non-zero linear partial difference 

operator with constant coefficients. There exists a function /(mi, . . . , m^) defined on Z" such 
that 

P/(mi,...,m„) = (5(mi,...,m„) . 

First Proof (Ehrenpreis-style) Consider the infinite-dimensional vector space, C[zi, . . . ,Zn,Zi^, ■ ■ ■ ,z~^], 
of all Laurent polynomials in zi, . . . , z„. Every function / on Z" uniquely defines a linear func- 
tional Tf defined on monomials by 

and extended by linearity. Conversely, any linear functional gives rise to a function on Z". Let 
P{zi, . . . , Zn) be the symbol of the operator V. We are looking for a linear functional T such that 
for every monomial z^^ ■ ■ ■ z^^" 

T{P{z,,. . . z^)zT' • • • C") = TsizT' • • • C") • 

By linearlity, for any Laurent polynomial a{zi, . . . , z„) 

T{P{zi,...,Zn)a{zi,... ,z„) =Ts{a{zi,...,Zn)) . 

So T is defined on the (vector) subspace P{zi, . . . , z„)C[zi, . . . , Zn,z^ , ■ ■ ■ , z~^] oiC[zi, . . . ,Zn,Zi , . . . , z~^] 
By elementary linear algebra, every linear functional on the former can be extended (in many 
ways!) to the latter. QED. 

Before embarking on the second proof we have to recall the notion of formal power series and 
more generally formal Laurent series. 

A formal power series in one variable z is any creature of the form 

oo 

Y^ aiz' . 

i=0 

More generally, a positive formal Laurent series is any creature of the form 

oo 

E 



OO 

OiZ" , 



where m is a (possibly negative) integer. On the other hand a negative formal Laurent series is 
any creature of the form 



E 

where r?2 is a (possibly positive) integer. 



oo 
i<.m, 



A bilateral formal Laurent series goes both "svays 

oo 

E 



OO 

^= — oo 



Note that the class of bilateral formal Laurent series is an abelian additive group, but one can't 
multiply there. On the other hand one can legally multiply two positive formal Laurent series by 
each other, and two negative formal Laurent series by each other, but don't mix them! Of course 
it is always legal to multiply any Laurent polynomial by any bilateral formal power series. But 
watch out for zero- divisors, e.g. 



(i-)E 



Any Laurent polynomial p{z) = UiZ^ + ■ ■ ■ ajZ^ of low-degree i and (high-)degree j in z (so a^ 7^ 0, 
Qj 7^ 0) has two multiplicative inverses. One in the ring of positive Laurent polynomials, and 
the other in the ring of negative Laurent polynomials. Simply write p{z) = z^aiPo{z) and get 
l/p{z) = z~'^{l/ai)po{z)~^ , and writing po{z) = 1 + qo{z), we form 



Po{z)-' = il + qoiz))-'=J2(-^y'io(' 



i=0 

and this makes perfect sense and converges in the ring of formal power series. Analogously one can 
form a multiplicative inverse in powers in z~^. 

It follows that every rational function P{z)/Q{z) in one variable, z, has (at least) two inverses, one 
pointing positively, one negatively. 

What about a rational function of several variables, P{zi, . . . , Zn)/Q{zi, . . . , 2;„)? Here we can 
form 2"'n! inverses. There are n\ ways to order the variables, and for each of these one can decide 
whether to do the postive-pointing inverse or the negative-pointing one. At each stage we get a 
formal one-sided formal Laurent series whose coefficients are rational functions of the remaining 
variables, and one just keeps going. 

Second Proof (Constructive): To every discrete function /(mi, . . . ,m„) associate the bilateral 
formal Laurent series 

^ /(mi,...,m„)zj" 

(mi,...,m„)GZ" 

We need to "solve" The equation 

p(zr\...,z-i)[ Y. /(mi,...,m„)zr---C"l =1 

\ (mi,...,m„)GZ" 



So "explicitly" 

{mi,...,m„)eZ" 

and we just described how to do it in 2"n! ways. 

The Maple package LEON 

This article is accompanied by a Maple package LEON. One of its numerous procedures is FS, that 
implements the above constructive proof. LEON can also compute polynomial bases to solutions 
of linear partial difference equations with constant coefficients, compute Hilhert Series for spaces 
of solutions of systems of linear differential equations, as well as find multiplicity varieties, in the 
style of Ehrenpreis, to 0-dimensional ones. 

Leon Ehrenpreis (1930-2010) A truly FUNDAMENTAL Mathematician (a Videotaped 
lecture) 

I strongly urge readers to watch my lecture, available in six parts from YouTube, and in two parts 
from Vimeo, see: 

http : //www . math . rutgers . edu/~zeilberg/mamariiii/inainarimhtnil/leon . html 

That page contains links to both versions, as well as numerous input and output files for the Maple 
package LEON. 
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